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Abstract —In this paper, we devise a scheduling algorithm for 
ordering transmission of synchrophasor data from the substation 
to the control center in as short a time frame as possible, 
within the realtime hierarchical communications infrastructure 
in the electric grid. The problem is cast in the framework of the 
classic job scheduling with precedence constraints. The optimiza¬ 
tion setup comprises the number of phasor measurement units 
(PMUs) to be installed on the grid, a weight associated with each 
PMU, processing time at the control center for the PMUs, and 
precedence constraints between the PMUs. The solution to the 
PMU placement problem yields the optimum number of PMUs 
to be installed on the grid, while the processing times are picked 
uniformly at random from a predefined set. The weight associated 
with each PMU and the precedence constraints are both assumed 
known. The scheduling problem is provabiy NP-hard, so we 
resort to approximation algorithms which provide solutions that 
are suboptimal yet possessing polynomial time complexity. A 
lower bound on the optimal schedule is derived using branch and 
bound techniques, and its performance evaluated using standard 
IEEE test bus systems. The scheduling policy is power grid¬ 
centric, since it takes into account the electrical properties of the 
network under consideration. 

Index Terms —Scheduling, NP-hard, approximation algo¬ 
rithms, hierarchical information flow, phasor data transmission. 

I. Introduction 

In a typical electric power grid, the communications infras¬ 
tructure is centered around communications between the indi¬ 
vidual substations and control centers. In the existing frame¬ 
work, this communications structure has disadvantages in that 
it offers slow automatic control by the control center and even 
slower manual control by system operators m. To overcome 
such drawbacks, installing GPS satellite-synchronized PMUs 
El across the grid which can record phasor data at very 
high frequencies are proposed. As per the year 2014, 1100 
PMUs had been installed across the U.S. Eastern Interconnect 
offering substantial coverage of the transmission system El. 
while China has been following large-scale implementation of 
PMUs for wide area monitoring ||4l. 

The high frequency of operation of PMUs results in accu¬ 
mulation of the voluminous phasor data at the substations for 
further processing and transmission. This is especially true 
with the proposition by the electric utilities to install large 
numbers of PMUs, resulting in overwhelming the existing 
communications network on the grid. The problem was high¬ 
lighted in a vision paper on the transmissions infrastructure 
for the future grid 0, where a hierarchical communications 
architecture for realtime wide area monitoring and control 
was examined. Basically, the following architecture was en¬ 
visioned; Phasor data at the substations form the first layer 
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Fig. 1. Time division multiplexing of PMU transmission. 

of information. Grid status routers form the interfacing layer 
between the substations and control centers, which form the 
higher level of information. Phasor data is transmitted from 
the substations to control centers via status routers over 
dedicated communication channels (see 0 Section I B]). This 
architecture is preferred over the more traditional redundant 
communications architecture (see 0 Fig. 3-3]), which presents 
problems related to management of phasor data especially with 
the increasing numbers of PMUs on the grid. 

In this paper, we base our study within the realtime hierar¬ 
chical communications infrastructure. Essentially, we concern 
ourselves with the following problem: Devise an algorithm to 
schedule transmission of phasor data (collected from a network 
of PMUs) from the substation to the control center in as short 
a time frame as possible, given the hierarchical information- 
flow architecture proposed in 0. An overview of the problem 
setup is presented in the following paragraph. 

We consider a power network with B buses and K branches. 
There are N < B PMUs in the grid, which could, for example, 
arise from optimal PMU placement algorithms 0 - |[T4l . Time 
is divided into frames with the duration of each frame equal 
to t units. A time frame is further divided into N slots each 
of duration ^ time units, as illustrated in Fig. [1] Within a 
time frame, the substation (which has gathered phasor data 
from the N PMUs) transmits to the control center via the grid 
status routers on dedicated channels of finite capacity. Given 
this setup, a transmission schedule for phasor data from the 
substation is proposed, so that the control center can use the 
ordered set of phasor data to more quickly and more reliably 
determine changes in the system state. Specifically, referring 
to Fig. [T] the fundamental question that we ask in this paper is 
the following: Would it be possible to minimize the processing 
time of phasor data from the N PMUs across the grid, i.e., to 
devise a scheduling policy that allows phasor data from the N 
PMUs to be transmitted from the substation in the specified 
order in as small a time frame t as possiblel 

The question posed in the previous paragraph is answered 
by casting the scheduling problem in the framework of single 
machine job scheduling with precedence constraints US). In 








simple words, we come up with a schedule to process phasor 
data from the optimal N PMUs in as small a time frame t as 
possible, under the precedence constraints between the PMUs; 
the precedence constraints are imposed to quantify the impor¬ 
tance or significance (in terms of electrical connectedness) of 
a PMU compared to other PMUs in the network. Since the 
problem is provably NP-hard, we address it via approximation 
algorithms, which provide solutions that are suboptimal yet 
possessing polynomial-time complexity; this component is 
the main contribution of the paper. The scheduling analysis 
presented herein is applicable without modifications to any 
configuration of PMU placement with precedence constraints 
of any form. 

We now discuss the role of precedence constraints in the 
scheduling problem. In the hierarchical information flow ar¬ 
chitecture, a network of PMUs transmit to a central monitoring 
station, where phasor data is archived for further processing. 
Employment of PMUs in backup protection schemes has been 
a popular practice ifThI - ifT^ . Consider a wide-area monitoring 
application with a fault detection and isolation scenario. It 
may be that faults on certain important transmission lines have 
much larger impact on the power grid than other transmission 
lines. In this case, it is important that the PMU data from the 
nodes corresponding to these transmission lines be reported 
to the central monitoring station earlier than data from other 
nodes. In another scenario, if phasor data from a particular 
PMU is found to be anomalous (transmission lines breaking, 
unexpected power demand trend, cyber attacks, unauthorized 
remote access to substation databases, etc.), then it might 
need to be assigned a higher precedence to be transmitted 
sooner to the control center (for timely monitoring and control) 
compared to phasor data from other PMUs in the network. In 
this case, the precedence constraints can change during the 
operation which is a more aggressive scenario than the one 
considered in this paper. The main point is that it is quite 
likely that application requirements will lead to precedence 
constraints on data transmission from various PMUs in the 
grid. In this paper, for analytical purposes, the precedence 
constraints are assumed known. 

In the first author and his colleagues devised a PMU 
scheduling scheme for transmission of phasor data from a 
network of N PMUs to the control center without a spe¬ 
cific communications architecture, as taken into account in 
the present paper. The objective there was to improve the 
performance of a fault detection scheme using scheduling, by 
accumulating phasor data from “all” the N PMUs to execute 
the fault detection procedure. This was clearly undesirable 
for realtime processing, since the scheme incurred a delay of 
{N — l)t/N time units per frame for each PMU which can 
trigger a hazardous situation even for moderate values of N 
and t. In this paper, we are motivated by realtime aspect of 
the hierarchical information-flow architecture proposed in Q. 
The problem considered in this paper is especially relevant 
when it is required to gather and analyze synchrophasor data 
for real-time decision making and wide area monitoring. 

Communications scheduling in power networks has been 
addressed from various angles. On the one hand, there are 
algorithms for power scheduling, while on the other hand. 


scheduling policies have been devised for multimedia data, 
user-access for smart power appliances, etc. For instance, a 
power scheduling scheme for the smart grid was proposed 
in fT\\ to improve the quality of experience (QoE); the QoE 
represents the customers’ degree of satisfaction. A multi-time 
scheduling scheme was presented in ll22ll in the framework 
of Markov decision processes for traditional and opportunistic 
energy users, under supply uncertainty due to variable and 
non-stationary wind generation, demand uncertainty owing to 
the stochastic behavior of a large number of opportunistic users 
and the coupling between sequential decisions across multiple 
timescales. A scheduling scheme for smart grid traffic (control 
commands, multimedia sensing data and meter readings) in a 
priority-based manner, aided by wireless sensor networks and 
cognitive radio technology, was presented in ll2^ . In ll24ll . a 
joint access and scheduling policy was proposed for in-home 
appliances (schedulable and critical) to coordinate the power 
usage to keep the total energy demand for the home below a 
target value. 

Contrasting to work in the existing literature, the focus of 
this paper is scheduling, or ordering, transmissions of phasor 
data from the substations to the control center on the power 
network by devising a feasible schedule in as small a time 
frame as possible. We are not concerned with power and grid- 
traffic scheduling that have been widely addressed by the smart 
grid community. To the best of our knowledge, this work and 
the one presented in ll20l are the first results on phasor data 
scheduling reported in the literature. 

The remainder of the paper is organized as follows. In 
Section HIl a mathematically precise problem statement is 
presented. Section |III] comprises the procedure to compute 
the weights associated with the PMUs, and the details of the 
approximation strategy to address the scheduling problem. The 
scheduling scheme is presented in Section |IV] An illustrative 
example along with the experimental results are provided in 
SectionjV] Conclusion and possible avenues for future research 
are provided in Section |VT] 

II. Problem statement 

Consider the set A/^ = N} of N PMUs installed on 

the power grid. The network of PMUs transmit phasor data to 
the substation where it will be archived for local processing 
and global operations (such as wide area monitoring, etc.) 
0. Phasor data from each PMU n G Af has a positive 
processing time, p„ > 0 and a nonnegative weight u>„ > 0; 
for purpose of analysis, the weights are assumed to be known. 
Let Cn denote the time at which phasor data from the n* 
PMU is processed completely by the control center in a 
feasible schedule. The scheduling scheme is required to satisfy 
precedence constraints, which are specified in the form of a 
directed acyclic graph Q = {Af,V) where {i,j) G V implies 
that the phasor data of PMUi must be completely processed 
by the control center before the substation can transmit data 
of PMUj. The main objective is to And a feasible schedule for 
the substation to transmit phasor data so as to minimize the 
weighted sum of the processing times, which is given by 

nGjV 


( 1 ) 


Preemption and control-center idle time are forbidden and 
hence phasor data all the PMUs should be processed in the 
interval where t = Preemption refers to a 

temporary interruption of processing, while idle time refers to 
the lapse between the end-time of processing phasor data of 
PMUi and the start-time of processing data of PMUj. 

In the classic setting, the aforementioned problem falls in 
the framework of single-machine job scheduling with prece¬ 
dence constraints ifTsll . The machine scheduling problem was 
shown to be strongly NP-hard ll25l . ll26l . and has received 
considerable attention both from a theoretical viewpoint EU 
- 1291 and for devising practical strategies l30l - l32l . Given 
the difficulty in devising efficient polynomial-time algorithms 
for optimally solving NP-hard problems, it is the usual practice 
to resort to polynomial-time suboptimal schemes, referred 
to as approximation algorithms 1^ . The approximation is 
optimal up to a small constant factor (for example, within 
2% of the optimal solution). An a-approximation algorithm 
runs in polynomial time and produces for every instance a 
feasible schedule of cost at most a times that of an optimal 
schedule. The value a is called the performance guarantee, or 
integrality gap, of the algorithm. With the increasing demand 
for installing a very large number of PMUs on the power 
grid, and given the societal importance of the problem, we 
will develop power grid-centric approximation algorithms to 
efficiently solve the scheduling problem posed in ([TJ. 

III. Solution methodology 

In this section, we elaborate on the approximation strategy 
to solve the scheduling problem posed in the optimization 
setup O. which is the subject topic of this paper. As men¬ 
tioned in the earlier sections, we assume knowledge of the 
weights Wn, n € J\f, and precedence constraints between the 
PMUs. 

We first present a 0-1 programming formulation of the 
optimization setup O- We define a variable Sum as follows; 

j 1, if PMU„ is processed before PMUm 

^nm ~ ^ ^ , . (^) 

10, otherwise. 

Let = 1 when the precedence constraints specify that 
PMU„ is a predecessor of PMUm and let 7„m = 0 otherwise. 
The processing of PMUm occurs at time 'Y^^Pn^nm + Pm- 
Therefore, ([T]i can be written as 


min PnSnmWm 

n m 

+ 

^ ^ Pm'^iTi 
m 

(3) 

subject to 5nm 

> 

^nm 

(4) 

^nm “1“ ^mn 

= 

1 

(5) 

^nm “1“ ^mk ^kn 

> 

1 

(6) 

^nm 

e 

{0,1} 

(7) 

^nn 

= 

0, 

(8) 

n = m = 

k = 

: 1,..., W, n ^ 

m, 

n ^ k, m ^ k. The constraint 0 specifies that Snm - 

= 1 


whenever PMU„ is a predecessor of PMUm, while (|5]l spec¬ 
ifies that PMU„ is to be scheduled either before or after 
PMUm. The matrix 5 = {5nm) is the adjacency matrix 


of a complete directed graph Qs- We introduce a constraint 
5mn + 5nk + ^km < 2 to ensure that the graph Qs contains 
no cycles. The constraint Q implies that the aforementioned 
constraint is equivalent to (|6]), which, therefore, ensures that 
Qs contains no cycles. When all the constraints are satisfied, 
Qs defines a complete ordering of the PMUs. The coefficient 
PnWm or 5nm in OJ denotes the cost of scheduling PMU„ 
before PMUm- The cost of scheduling PMU„ before PMUm 
is denoted in terms of the cost matrix C = {cnm), where 

_ j Pn^m: if ^mn — 0; 

Cnm — S -r- 1 

I^CX), if '~frnn — I- 

Whenever the precedence constraints specify that PMUm is a 
predecessor of PMU„, we have c„m = oo. The setup (O can 
now be written as 


min EE Cnm^nm H” E Pm^mi (lO) 

n m m 

subject to constraints specified by (|5]l, (|6]l, (|2]l and ®. Each of 
the constraints in (|5]l is incorporated into a Lagrangian function 
with multipliers anm as follows 


= EE( 


Cnm ^nm ^mn J 


)s„ 


EE ^ ^ Pm'^m ? 

n m^n m 


(ii) 


which is to be minimized subject to Q, ([Hi, ([tIi and We 
let 


^nm — ^mn — ^ min^C 


nm? ^mn 




( 12 ) 


so that the multipliers provide as large a contribution as pos¬ 
sible to the lower bound. With this, we now define a reduced 
cost matrix = {c^nm), where C^nm = C„m - - Umn- 

In the presence of cycles, constraints © and (|6l) can be 
used to derive cycle elimination constraints involving q edges 
of the form 


9 

i=l 



(13) 


where ni = Uq+i and ni,...,nq correspond to q differ¬ 
ent PMUs. Suppose that a Lagrangian relaxation of r — 
I of the constraints (O are performed using multipliers 
giving the following Lagrangian function 


^7-1) ^ + 

n m n m^n 

r — 1 

k—l m 

The constraint (fOl l is introduced into (fT4li to get 

L(r) = ^ (15) 

= , min 1. 

Choosing /3L) as large as possible retains the nonnegativity of 
the coefficients of the variables. The update equation for the 


reduced cost matrix is = C^'' where the 

matrix = (bnm) is defined as follows: 




1, if n = rii and m = rii+i, 
0, otherwise, 


(16) 


where i = 1,..., q. Thus, we can write in the same form 
as L(’-i) as given by (Ell. Using L^‘^\ a lower bound on the 
optimal solution of ([T]i can be obtained. 

Theorem 3.1: The following is a lower bound for the sum 
of weighted completion times: 


LbM=^ ^ + (17) 

n m^n 2=1 m 

Proof: From the aforementioned development, we see that 
Cnm > 0, n,m = 1,..., N and n m. Therefore, given 
by (fTSl l can be minimized by setting 6nm = 0 whenever Cnm > 
0 which implies Cnm^nm = 0. This proves Theorem 13. II ■ 

From (O, it is clear that the lower bound comprises 
only cycle elimination constraints. Furthermore, increasing the 
lower bound amounts to introducing more such constraints 
which is equivalent to increasing the multipliers. In other 
words, to obtain the tightest lower bound, as many cycle 
elimination constraints should be introduced. The following 
theorem states that when no more such constraints can be 
found, the reduced cost matrix can be used to find a 
feasible schedule of the PMUs. 

Theorem 3.2: There exist variables satisfying the con¬ 
straints (Hli - (Isll such that CnmSnm = 0 if and only 

if no additional constraint with a positive multiplier can be 
introduced into the Lagrangian function given by (flST l. 

Proof: Consider a directed graph Q, where a vertex 
denotes a PMU, and if Cnm > 0 then there exists an edge 
between the vertex pair {n,m). The variables satisfying dUl 
- (HI) and = 0 define an ordering of the 

PMUs for which the order graph Qs contains no edge of Q. 
In other words, the graph obtained by reversing all edges 
of Os contains all edges of Q. This implies Q contains a 
partial order without cycles; therefore, no cycle elimination 
constraints with a positive multiplier exists. On the other hand, 
if no cycles exist in 0, then there exists an ordering consistent 
with the partial ordering in Q. The reverse ordering provides 
the values of the variables which satisfy (0]) - ([8]) such that 
CnmSnm = 0. TWs proves Theorem |32] ■ 

If a schedule of PMUs can be found from by defining 
values of the variables that satisfy the constraints (0]) - ([8]), then 
the sum of the weighted completion times for this schedule is 
an upper bound computed as follows: 

Ub(0 = LB(’-) + ^^ - l\ , (18) 

i—l \ n m / 

where the matrix B^), i = 1,..., r is defined by (fTbl l. 

Corollary 3.3: LB^’’^ = UB^*"^ if and only if all constraints 
in the Lagrangian function are satisfied with equality. 

From Corollary 13.31 we see that the lower bound can be 
made tighter by using the constraints that are satisfied as 
inequalities. Let the PMUs be renumbered so that the schedule 


obtained from is (1,...,A^), and therefore, is a 

lower triangular matrix. Now, suppose that a constraint with 
multiplier /3 satisfying the strict inequality has been found. If 
this /3 is removed from the Lagrangian function, the entries of 
the matrix increase by a factor of /3. Thus, for n < m, 
each 6nm produces an element Cnm = P placed above the 
leading diagonal in C^'^\ New cycle elimination constraints 
can be introduced into the Lagrangian function by using such 
5nm for which Cnm > 0 along with Sij, n < j < i < m for 
which cl 3 > 0. If new cycle elimination constraints having 
multipliers which sum to 0 < /3* < /I are introduced into 
the Lagrangian function for each of 77 > 2 variables Snm 
of the original constraint, then the original chosen constraint 
is reintroduced into the Lagrangian function with a reduced 
multiplier /3 —/3*, and the lower bound increases by (p —1)/3*. 
By choosing rj and /3* as large as possible, we can provide the 
maximum increment to the lower bound. Note that, the new 
Lagrangian function necessitates finding a new schedule such 
that is a lower triangular matrix. 

To find new constraints involving the variable S^m {n < w) 
mentioned in the aforementioned paragraph, we present the 
following simple procedure: Consider a network of buses 
n,... ,m, where m and n denote the source and sink, respec¬ 
tively. An edge exists between buses i and j with capacity 
clj'^ whenever n < j < i < m. A capacity of P is placed 
on the source bus. The problem of finding the maximum 
flow from the source to sink buses is equivalent to that 
of generating new constraints. The maximum flow can be 
decomposed into different flows along paths from source to 
sink buses. The variables corresponding to the edges in a path 
from source to sink can be combined with Snm to produce 
cycle elimination constraints; the associated multiplier is given 
by the flow along that path. When is reduced, the entries 
remain nonnegative, since the flow along any edge cannot 
exceed its capacity. Also, the lower bound is the best possible 
for the constraints obtained from the network flow problem. 
This is because, all constraints contain at least one variable 
corresponding to an edge which cuts across the minimum cut 
set, and the fact that the maximum flow is equal to the value 
of the minimum cut set. 

A. Method to achieve the lower bound 

The first step is to compute given by (fTTT i which 
requires 0{N‘^) computations. The second step is to decide 
the cycle removing constraints and the order in which they are 
to be introduced into the Lagrangian function. 

For the bound to be tight, from Corollary 13. 3l we know that 
any constraint that is added should be satisfied with equality, 
and is a suggestive rule indicating which constraint should be 
used. Given this condition, for small values of q inequality 
(fT3l) is satisfied with equality. Also, a schedule generated by a 
heuristic procedure indicates whether the constraint is satisfied 
as an equality. With this setup, we add all cycle elimination 
constraints with three edges in 0{N^) steps that are satisfied 
as equalities when the values of the variables correspond to 
heuristic schedule. This process is repeated for constraints with 
four edges in 0{N'^) steps. 


Lastly, a technique is needed to remove remaining cycles 
and to enable a schedule to be found from the reduced cost 
matrix It is easy to see that any PMU n with Cnm = 0 
for all PMUs m can be placed in the first unfilled position in 
the schedule, and row n and column n can be deleted from 
the reduced cost matrix. If there is a choice of PMUs that can 
be scheduled, the one that is scheduled first in the heuristic 
method is chosen. The process if repeated until no PMU can 
be scheduled in the first unfilled position of the schedule. In 
the event that a complete schedule has not been generated, 
it is possible to schedule PMUs in the final position of the 
schedule in a similar manner. 

From Theorem 13.21 it follows that when unscheduled PMUs 
remain, at least one additional constraint can be included 
in the Lagrangian function to tighten the lower bound. We 
pick a PMU„ such that the number m of unscheduled PMUs 

fr') 

with Cnm > 0 is as small as possible. A search procedure 
involving 0{N'^) steps can be invoked if a cycle elimination 
constraint with PMU„ exists. In the absence of a cycle, another 
unscheduled PMU is similarly examined, and the process is 
repeated till all the PMUs are scheduled. 

During implementation of the lower bound, we saw that 
0{N'^) steps are required for the computation of the lower 
bound. For the elimination of the cycles, the network flow 
problem, whose solution requires 0{N^) steps, is solved 
for each variable taking the value one in each constraint 
satisfied as a strict inequality. There are 0{N‘^) constraints 
each comprising 0{N) variables. Therefore, we need 0{N^) 
steps if we are not forced to generate a new schedule during 
the network flow problem. Thus, it is clear that the number 
of steps required by the network flow is pseudopolynomial, 
provided a new schedule is generated so that the lower bound 
keeps increasing by one unit. 

IV. Scheduling scheme 

In this section, we describe the heuristic method used before 
the branch and bound algorithm, along with the branching rule 
and the search procedure. 

We employ a simple heuristic method to obtain a schedule. 
Heuristics are used to indicate which constraints are to be 
used in the lower bound, and are not concerned with providing 
upper bounds to the optimization setup. Basically, we select a 
schedule at random that satisfies all the precedence constraints. 
To improve the heuristic, a group of PMUs is selected which 
consists of PMUs that are scheduled in consecutive positions 
in the given schedule such that there is a precedence constraint 
between each adjacent pair of PMUs. The PMUs of this group 
are removed from their current positions in the schedule and 
inserted elsewhere such that the new schedule satisfies all the 
precedence constraints in addition to having a smaller sum of 
weighted completion times than the previous sequence. The 
process is repeated until no further improvements are possible. 

The central idea behind the branching rule is to reduce 
the difference between the lower and upper bounds computed 
at the bus from which we decide to branch. We search for 
constraints which satisfy the inequalities and the one which 
has the largest multiplier is selected. The variable 5nm which 
appears in this constraint is selected so that 5nm takes the value 


one, and such that PMU„ is not a predecessor or a successor 
of PMUm- The search tree with two branches are formed, such 
that in one branch PMUm is constrained to be scheduled be¬ 
fore PMU„ and in the other branch PMU„ is constrained to be 
scheduled before PMUm. The precedence constraints between 
other pairs of PMUs are directly implied by transitivity. The 
variable 5ij is removed from the Lagrangian without altering 
the lower bound if PMU^ is scheduled before PMUj. In the 
first branch of the search tree, where PMU„ is constrained 
to be scheduled after PMUm, new schedules will emerge 
after the addition of more cycle elimination constraints. In the 
second branch, where PMU^ is constrained to be scheduled 
before PMUm, the constraint from which the variable 5nm is 
originally found will be removed. 

The goal of the search strategy is to specify from which 
bus to branch. A simple active node search is employed 
which uses backtracking to compute the reduced cost matrix 
from and the list of constraints. Doing so, we avoid 

storing the reduced cost matrix at node of the search tree. 
In the following section, we perform computer simulations to 
validate the theoretical findings of this paper. 

V. Illustrative Example 

The algorithm was tested on different grids. For each 
grid, the optimal number of PMUs was obtained by solving 
the PMU placement problem, with the objective to achieve 
complete network observability ignoring zero injection mea¬ 
surements. For phasor data from each PMU, an integer pro¬ 
cessing time (in ms) from the uniform distribution [1, 50] is 
chosen. The algorithm was implemented using MATLAB® on 
a Windows 7 PC, with a dual-core central processor, and 4.00 
GB random access memory. 

In the theoretical analysis carried out in the earlier sections 
of the paper, we had assumed knowledge of the weights and 
the precedence constraints. For the purpose of the experimental 
study presented in this section, we will utilize a simple 
procedures to derive the weight Wn, n = 1 ,..., of each 
PMU and the precedence constraints between the PMUs. This 
is being done for this illustrative example. As noted earlier, 
we expect that the precedence constraints will come naturally 
from the power grid application requirements. 

The SVD of the bus admittance matrix is a useful measure 
of the electrical connectedness between various components 
in the power network. A thorough exposition on this topic 
was first reported in El, where the authors introduce the 
term “electrical centrality measure” to quantify the degree of 
connectedness between buses in the grid. It was analytically 
demonstrated in B1 that typical grid admittance matrices have 
singular values and vectors with only a small number of strong 
components, and performing SVD on the admittance matrix 
leads to no loss of connectivity-information. 

The weights Wn, n G AA are obtained as follows: 

1) We obtain the optimal number {N) of PMUs by solving 
the PMU placement problem. 

2) In the B x B bus admittance matrix, we choose the rows 
and columns which correspond to the bus numbers where 
PMUs are installed, leading to an W x sub-matrix. 


3) We perform the SVD of the N x N sub-matrix to obtain 
the singular values and singular vectors. The x 1 left and 
right singular vectors are denoted and Vn, respectively, 
while the singular values are denoted an- 

4) We compute the Euclidian norm of the vectors and 

define the weight associated with the PMUs as follows: 


A 

Wn = 


N 




(19) 


Note that, the index of each entry of the vector 
corresponds to a bus location where a PMU is installed. 


The precedence constraints are obtained as follows: 

1) In the vector aiUi, i.e., the first column of the N xN sub¬ 
matrix obtained as above, the PMU placed on the entry with 
the highest magnitude is given the highest precedence. Ui 
is the eigenvector corresponding to the largest eigenvalue. 

2) The procedure in Step [U is repeated for the remaining 

vectors n = 2,...,N, where the m„s are picked 

in the decreasing order of the corresponding eigenvalues. 

However, in the aforementioned procedure there is a possibil¬ 
ity of conflict. For example, consider two vectors aiUi and 
a^u^. Suppose the entry having the largest magnitude in vector 
aiUi is the same as the entry having the largest magnitude 
in vector a^u^. Then, the procedure picks the same PMU in 
both the vectors ctiMi and a^u^. To resolve this conflict, we 
propose the following modification: Pick the entry in the vector 
(73 W 3 having the second largest magnitude. If this entry is not 
the same as the one in vector (T 2 M 2 , then the PMU placed on 
that entry is given a lower precedence than the PMU obtained 
in the vector a 2 U 2 - This simple procedure is implemented for 
all the vectors anUn- 

We consider the IEEE 14-bus system to illustrate the 
procedure adopted for establishing the precedence constraints 
between the PMUs. We obtain the optimal number PMUs 
employing the topology of the grid cni, im. We consider 
a single time frame, and implement the following steps: 


Bus number 
2 

6 

7 

9 



12.8874 


22.7655 


43.1198 


5.7771 


I 


14.7761 


4.7641 


8.6612 


1.2276 


(TiXli 0‘2U2 Cr3U3 (T4U4 


Time frame t 

Fig. 2. Precedence between PMUs for the IEEE 14-bus network. The values 
of the elements of the vectors anUn are also indicated. 


complete network observability ignoring zero injection 
measurements (see n, Qo), ini). 

(2) In the 14 X 14 bus admittance matrix, pick the rows and 
columns numbered 2, 6, 7 and 9, thereby yielding a 4 x 
4 sub-matrix. 

(3) Perform the SVD of the 4x4 sub-matrix to obtain the 4 x 

1 right and left singular vectors and i;„, respectively, 
and the singular values cr„, n = 1, ..., 4. Compute the 
magnitude of the elements of the vectors The index 

of each entry of the vector corresponds to a bus 

location where a PMU is installed. The four vectors 

n = 1,..., 4 are depicted in Fig.|2] where a column denotes 
a vector, while a box in each column denotes an entry of 
the vector. The magnitude of the elements of vectors 
are also indicated. The number of boxes in each column 
equals the number N of PMUs installed on the bus system. 

(4) In the vector aiUi (the first column in Fig. |2]i, the entry 
having the largest magnitude appears in the last row - 
marked in blue. Thus, the PMU placed on bus numbered 
9 is given the highest precedence. 

(5) In the vector (T 2 M 2 (the second column), the entry having 
the largest magnitude appears in the last row, similar 
to that in the vector ctiMi, again allocating the second 
highest precedence to the PMU placed on bus numbered 
9 (conflict). However, this conflict is resolved by giving 
the second highest precedence to the PMU placed on the 
bus numbered 6, which has the second largest magnitude 
in the vector a 2 U 2 - In other words, PMUg > PMUg. 

(6) Continuing in this fashion, and employing the conflict- 
resolution strategy, the PMUs installed on buses 2, 6, 7 
and 9 are given the precedence in that order, i.e., PMUg > 
PMUs > PMUt > PMU 2 . 

We thus have precedence constraints of the following form 
between the PMUs: For n > m, PMU„ > PMUm, i-e., PMU„ 
is designated to be transmitted before PMUm- 

We consider the IEEE-14, IEEE-30, IEEE-39, IEEE-57, 
IEEE-118, IEEE-300 and power flow data from two Polish 
electric grid systems. The bus and line data for the bus 
networks were obtained using MATPOWER® ll35l . In Table H] 
we show the different bus networks considered in this paper 
along with the number of branches and the optimal number 
(N) of PMUs to be installed on each of the grids for complete 
network observability ignoring zero injection measurements. 


Bus 

network 

Number of 
branches 

Optimal (A^) 
PMUs 

IEEE-14 

20 

4 

IEEE-30 

41 

10 

IEEE-39 

46 

13 

IEEE-57 

80 

17 

IEEE-118 

186 

32 

IEEE-300 

411 

156 

Polish-2737 

3506 

971 

Polish-3375 

4161 

1384 


TABLE I 

Bus NETWORKS CONSIDERED IN THIS PAPER AND THE OPTIMAL N PMUS 
OBTAINED BY SOLVING THE PMU PLACEMENT PROBLEM (??). 


(1) For the 14-bus network, an optimum of W = 4 PMUs 
is to be placed on buses numbered 2, 6 , 7 and 9 for 


We compare the performance of our scheduling scheme with 
that of a simple greedy algorithm. We begin by considering 
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a directed acyclic graph based on the precedence constraints 
which are assumed known. Then, we pick that PMU from the 
directed acyclic graph which has no parent to transmit, and 
then delete that PMU from the graph. If there are more than 
one PMU, we pick the one with the largest weight. 

In Fig. [3 we show the plot of average computation time as a 
function of the number of PMUs installed on the grid for both 
the greedy algorithm and for the scheduling scheme devised in 
this paper. We see that, for less than 200 PMUs, our scheduling 
scheme has good performance with regard to computation 
time. However, when more number of PMUs are installed, the 
average computation time increases significantly. The reason 
for an increase in the average time for larger number of PMUs 
is basically computational in nature. Also, the performance 
of the scheduling scheme is comparable to that of the greedy 
algorithm for fewer number of PMUs. As the number of PMUs 
increases, the greedy algorithm behaves very poorly, which 
corroborates intuition; greedy algorithms approximate a global 
optimal solution in reasonable time for fewer instances of the 
problem, however, as the number of instances increases, the 
performance degrades. A similar behavior is observed in Fig.|4] 
when we plot the average number of nodes required in the 
search tree as a function of N. 



Fig. 3. Average computation time versus the optimal number [N) of PMUs, 
using different computational resources. 



Max. computation 
time (sec) 

Avg. computation 
time (sec) 

Bus 

network 

Optimal (A^) 
PMUs 

GA 

AOA 

GA 

AOA 

IEEE-14 

4 

0.182 

0.312 

0.061 

0.091 

IEEE-30 

10 

0.371 

0.598 

0.112 

0.132 

IEEE-39 

13 

0.612 

0.861 

0.158 

0.198 

IEEE-57 

17 

1.125 

2.761 

2.118 

1.41 

IEEE-118 

32 

18.112 

5.871 

3.547 

3.87 

IEEE-300 

156 

23.129 

11.651 

7.242 

7.182 

Polish-2737 

971 

- 

- 

343.17 

277.12 

Polish-3375 

1384 

- 

- 

591.71 

391.65 


TABLE II 

Maximum and average computation times. GA - greedy 

ALGORITHM, AOA - APPROXIMATELY OPTIMAL ALGORITHM DEVELOPED 
IN THIS PAPER. 



Fig. 4. Average number of nodes in the search tree versus the optimal number 
(N) of PMUs. 

In Table [HI we present the results of maximum and average 
computation times for the greedy algorithm and the proposed 
scheduling scheme. Here again, we see that, as the number of 
PMUs increases, the maximum computation time to obtain a 
feasible schedule increases, owing to the size of the problem. 
It is also of interest to note that, for the last two entries in the 
Table [III one cannot obtain the maximum computation time, 
since the average computation time and the average number 
of nodes that are calculated and plotted in Fig. [3 and Fig. |4] 
represent lower bounds on the respective averages because 
algorithm does not converge in polynomial time. 

To demonstrate the impact of computational resources on 
the performance of scheduling, we also tested our algorithm 
on quad-core and six-core processors with 6GB of random 
access memory. The resulting average computation time for 
varying N is shown in Fig. [3 where we notice reduction 
in computation time as a function of the number of pro¬ 
cessors. We note that, our experiments are performed under 
laboratory conditions and is meant for illustrative purposes. 
In practice, the command center is endowed with far more 
powerful computational tools, which are indeed necessary for 
realtime operations. Using such tools, the performance of our 
scheduling scheme can be significantly improved, and can be 
fine-tuned to satisfy the stringent requirements on latency for 
transmission applications on the grid (see, for example, ESI). 
In order to avoid obscuring the main topic of the paper, we 
have not investigated computational implementation in greater 
detail. This work is a first step in the direction of devising 
scheduling algorithms for phasor data transmission, with scope 
for improvement in both theory and practice. 

VI. Conclusion and future work 

We have devised a scheduling algorithm for ordering trans¬ 
mission of phasor data from the substations to the control 
center in as short a time frame as possible on the hier¬ 
archical information architecture in the electric power grid. 
The problem was cast in the framework of classic single 
machine job scheduling with precedence constraints. The 
optimal number of PMUs was obtained by solving the PMU 













































placement problem with the objective of complete network ob¬ 
servability ignoring zero injection measurements. The weights 
associated with the PMUs and the precedence constraints 
between the PMUs were assumed known. The processing 
time of the PMUs were chosen uniformly at random from 
a pre-specified set of times. The problem is provably NP- 
hard, so approximation algorithms are typically used to obtain 
polynomial time complexity albeit being suboptimal. Using 
branching and bounding techniques, we arrive at a lower bound 
on the feasible schedule. The average computation time and 
the average number of nodes in the search tree was plotted 
as a function of the number of PMUs installed on the grid, 
and were compared with the performance of a simple greedy 
algorithm. As expected, for fewer PMUs on the grid, the 
performance of scheduling is favorable, however, for more 
than 200 PMUs, the computation time increases at an alarming 
rate. Future work would involve (i) devising more efficient 
cycle elimination techniques to reduce the gap between the 
lower and upper bounds and (ii) problem formulation for 
distributed scheduling to facilitate distributed functionality of 
the existing control centers ||5] Section I D]. 
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